
THEORY OF TWO-POINT 
CORRELATIONS OF JET NOISE 

Herbert S. Ribner 

Langley Research Center 
Hampton , Va. 23665 


« 





I' 

h. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • DECEMBER 1976 


TECH LIBRARY KAFB, NM 



ERRATA 




\|Page 

\]page 

£age 


^Page 

jPage 

^Page 


JPage 

Jp 


age 


'Wge 

''Wge 

Page 


NASA Technical Note D-8330 

THEORY OF TWO- POINT CORRELATIONS OF JET NOISE 


Herbert S. Ribner 
December 1976 


12, equation (15): 

14, equation (22): 
20, equation (38): 


Replace (u*) 2 (u£) 2 with (u*) (u^f 
Replace dflrf with d^if d^Y. 

The first line of the right side of the equation should read 


nanb^anb^ 2 


4,2c 0 4 x 2 JJ f 

20, line preceding equation (40): Replace U^(y) with uj^(y). 

21, equation (42): Replace the factor g 2 g 3 with th^ factor p^2 g2&3* 

22, equation (43): Replace the factor with the factor — 

2 1/2 tt 2 1 /2 7r 

g 9 go p 2 g 9 go 

22, equation (44): Replace the factor - with the factor ° 


47T 


4ff 


26, line 3: Replace the word "larger" with the word "smaller" so that the 

sentence reads 

This yields manyfold smaller values of ff. 

27, line 12: Replace equations (16) with equation (32). 

29, line 10: Replace figures 7 and 10 with figure 7 and other figures. 

30: Move (and amplitude) on line 3 from bottom of page to end of sentence 
so that the sentence reads 

The sound emission from each small volume is seen to consist of a 
rather directional quadrupole pattern whose lobes fluctuate randomly in 
direction (and amplitude). 


ISSUED MAY 1977 



TECH LIBRARY KAFB, NM 


□ 13M 


23 


1. Report No. 2. Government Accession No. 

NASA TN D-8330 

3. Recipient's Catalog No. 

4. Title and Subtitle 

THEORY OF TWO -POINT CORRELATIONS OF JET NOISE 

5. Report Date 

December 1976 

6. Performing Organization Code 

7. Author(s) 

Herbert S. Ribner 

8. Performing Organization Report No. 

L-11058 

10. Work Unit No. 

505-03-11-02 

9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, VA 23665 

11. Contract or Grant No. 

13. Type of Report and Period Covered 
Technical Note 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 

14. Sponsoring Agency Code 

|l5. Supplementary Notes 


Herbert S. Ribner: Staff Scientist, now back at University of Toronto, Institute for 

Aerospace Studies. 


16. Abstract 

A large body of careful experimental measurements of two-point correlations of far- 
field jet noise has been carried out and was briefly reported recently by Lucio Maestrello in 
NASA TM X- 72 83 5. The results were not altogether as expected and motivated the present 
effort to bring theory to bear. The model of jet-noise generation is an approximate version 
of an earlier work of Ribner, based on the foundations of Lighthill. The model incorporates 
isotropic turbulence superimposed on a specified mean shear flow, with assumed space -time 
velocity correlations, but with source convection neglected. The particular vehicle is the 
Proudman format, and the previous work (mean -square pressure) is extended to display the 
two -point space -time correlations of pressure. 

The shape of polar plots of correlation is found to derive from two main factors; 

(1) the noncompactness of the source region, which allows differences in travel times to the 
two microphones - the dominant effect - and (2) the directivities of the constituent quadru- 
poles - a weak effect. The noncompactness effect causes the directional lobes in a polar 
plot to have pointed tips (cusps) and to be especially narrow in the plane of the jet axis. In 
these respects, and in the quantitative shapes of the normalized correlation curves, results 
of the theory show generally good agreement with Maestrello's experimental measurements. 


17. Key Words {Suggested by Author(s)) 

Jet noise 

Acoustics 

Aeroacoustics 


18. Distribution Statement 

Unclassified - Unlimited 


Correlations 


Subject Category 71 

19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 


58 

$4.25 


For sale by the National Technical Information Service, Springfield, Virginia 22161 




i i ii 1 1 ii i i i i i ii ii i mi i i 

CONTENTS 

SUMMARY 1 

INTRODUCTION 1 

SYMBOLS 3 

OVERVIEW OF THEORY 7 

ANALYTICAL DEVELOPMENT 8 

Governing Equations 8 

Self-Noise 12 

Integral Over Noncompact Jet 16 

Shear Noise 19 

Summarized Predictions 22 

DISCUSSION 24 

Comparison With Experiments of Maestrello 24 

Further Considerations 28 

CONCLUSIONS 30 

APPENDIX A - QUADRUPOLES WITH FLUCTUATING ORIENTATION: THE 

BEAM PATTERN CONCEPT 32 

APPENDIX B - EVALUATION OF NONCOMPACTNESS FACTORS FOR CASE I: 

0 a = 0 b =0 > ^a = 0 35 

APPENDIX C - EVALUATION OF NONCOMPACTNESS FACTORS FOR CASE n: 

<?b = < ?a 37 

APPENDIX D - REDUCTION OF THE SHEAR-NOISE INTEGRAL, 

EQUATION (42) 39 

REFERENCES 41 

FIGURES 43 


iii 



THEORY OF TWO-POINT CORRELATIONS OF JET NOISE 


Herbert S. Ribner* 

Langley Research Center 

SUMMARY 

A large body of careful experimental measurements of two -point correlations of 
far -field jet noise has been carried out and was briefly reported recently by Lucio 
Maestrello in NASA TM X-72835. The results were not altogether as expected and 
motivated the present effort to bring theory to bear. The model of jet-noise generation 
is an approximate version of an earlier work of Ribner, based on the foundations of 
Lighthill. The model incorporates isotropic turbulence superimposed on a specified 
mean shear flow, with assumed space-time velocity correlations, but with source convec- 
tion neglected. The particular vehicle is the Proudman format, and the previous work 
(mean-square pressure) is extended to display the two-point space-time correlations of 
pressure. 

The shape of polar plots of correlation is found to derive from two main factors: 

(1) the noncompactness of the source region, which allows differences in travel times to 
the two microphones - the dominant effect - and (2) the directivities of the constituent 
quadrupoles - a weak effect. The noncompactness effect causes the directional lobes in 
a polar plot to have pointed tips (cusps) and to be especially narrow in the plane of the 
jet axis. In these respects, and in the quantitative shapes of the normalized correlation 
curves, results of the theory show generally good agreement with Maestrello' s experi- 
mental measurements. 


INTRODUCTION 

A quarter of a century of jet-noise research has focused almost exclusively on the 
response of a single microphone. Until now there have been no studies of the joint 
response of two microphones in the far field, that is, of two-point correlations of sound 
pressure. This is understandable because the difficulty of theoretical prediction has been 
so formidable that one leans toward the simpler case. 

Recently, however, rather different theories (e.g., those of Mani (refs. 1 and 2) 
and Ribner (refs. 3 and 4)) have demonstrated considerable success in predicting 
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single -microphone directional response (refs. 2, 5, and 6). This points up a certain 
insensitivity of the single-microphone response (mean-square pressure) to details of the 
theoretical model. The two-microphone correlations, which provide much more informa- 
tion, would be expected to be less ambiguous in their dependence on details of the model. 
Such cross-correlations could conceivably provide a sensitive test of validity. 

With this motivation, a large body of careful experimental measurements of two- 
point jet-noise correlations has been carried out by Maestrello and is reported in a 
recent paper (ref. 7) along with theoretical considerations. The results were not 
altogether as expected and, in accord with the above notions, motivated the present 
attempt at a theoretical explanation. 

The model of jet-noise generation is an approximate form of that of Ribner (refs. 3 
and 4), which is a development of the Lighthill (refs. 8 and 9) formalism. It postulates 
isotropic turbulence superimposed on a specified mean shear flow. The Proudman 
(ref. 10) formalism is used, wherein the momentum flux in the direction of the observer 
governs the combined emission of the Lighthill quadrupoles in that direction. The space- 
time velocity correlations associated with this momentum flux are described by a simple 
mathematical model. The development of reference 4 (mean-square pressure) is extended 
to provide the two-point correlations of pressure. A series of simplifying assumptions 
and restrictions, based on physical considerations, is applied to render the sixfold inte- 
gration tractable. 

An important restriction is the complete omission of all source convection and 
refraction effects. The convection effects are expected to be predominantly amplitude 
changes at points of the sound field, with much less change in the pattern coherence 
between two points. Hence, normalizing the two-point pressure correlation by the product 
of the rms pressures at the two points should largely cancel the contribution of convection. 
The refraction effects can be anticipated with less assurance; nevertheless, the normali- 
zation should partially offset them. 

The single-microphone case likewise calls for a sixfold integration. However, the 
major features of the directivity can be inferred after only a threefold integration. This 
yields the sound emission from a unit volume of a jet whose directional properties are 
fairly representative of those of the entire jet in a given narrow frequency band. That is 
not so, however, for the two-microphone correlation; hence the much more demanding 
sixfold integration is required. 

The details of the derivation of the two-point pressure correlations are developed 
in the body of the paper and in appendixes A to D. The results are relatively simple 
closed-form expressions. The theoretical predictions are compared with the measured 
data of Maestrello (ref. 7) for a number of cases. 
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SYMBOLS 



f f 

Hi) 


g ( T i) 


Sl’®2’®3 

I,J 

L 


e. 


N Str 


P 

P 


a - y: 

space factor in correlation u?u° , specified in equations (21) 
ambient speed of sound 
jet diameter 

function defined in equation (27) 

frequency characterizing the local turbulence, treated as a "typical" frequency 

2 2 / T 2 

= e~ v 71 / L 

q Vj " CtXp I T ■t I 

time delay factor in correlation u^Uj , specified as e 1 A| 

constants 

noncompactness integrals appearing in equations (46) and appendix B 

effective scale of turbulence 

effective length of source region 

jet Mach number, Uy/c Q 

Strouhal number, f^I^/Uj 

unit vector along x 

ith component of n in basic reference frame (fig. 1) 

ith component of n in special reference frame (eqs. (17)) 

pressure 

far -field pressure perturbation; also point of pressure measurement 
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q se = 4 < w f> ( R / C o) | sin 0 b sin (V 2 )| 

q sh = 2 ( w f)( R / c o)| sin 0 b sin (V 2 )| 

R radius of jet; also correlation function 

r radial distance from jet axis 

T-. = o v.v. + r.. + IP - c ^pi 6-- 

1J F o 1 J T 1] \ o 7 lj 

t time of signal reception 

t time of signal emission 

U local mean flow velocity 

Uj jet nozzle velocity 

U x component of U along x 

n local turbulence velocity 

Uj component of n 

u x component of if along x 

v local resultant velocity, u + U 

Vj component of v 

v x component of v along x 

if position vector of field point 

Xj component of if 

Y position vector of centroid of correlation volume, or Tt eddy Tf 
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Y i 


component of Y 


P 

y 

A 


6 u 


position vector of source point 

_ 2 x Maximum shear -noise intensity ^ gee (45)) 

Self-noise intensity 

azimuth angle in cylindrical polar coordinates 

time delay between emissions for simultaneous reception at p a and p^ 

£ < w f> i 

= — - — cos e b - cos e a 


Kronecker delta, 6^ = (i ^ jj 


dummy variable 


V 

8 


separation of source points, y a - y^ 


component of if 

polar angle, x^ = x cos 6 

wavelength 


P 

P 

Pn 


azimuth angle in special reference frame (see sketch (a) following eqs. (17)) 


density 


ambient density 


constant (eq. (35)) 

a b 

time delay between signal receptions at x and x (specified) 


time delay between signal emissions, r + A 
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<p 


viscous compressive stress tensor 
azimuth angle, X 2 = x cos 9 cos cp 


\f/ angle between x a and x b ^Case I: cos xjy = cos 2 0 + sin 2 0 cos cp^ 

Case II: ^ | - ® a | ) 

O) radian frequency, 2rrf 

= 2 rfj 

Subscripts or superscripts: 

—•* 3 , 

a associated with sound emission to point x 

b 

b associated with sound emission to point x 


Subscripts: 

i ith vector component (e.g., i = 1 denotes component along xj) 

i,j tensor indices 

n radial direction (normal to shear layer) 

se self-noise 

sh shear noise 

Brackets: 


CD 

<> 


evaluated at retarded time 
average over mixing region 
A bar over a quantity denotes a time average. 
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OVERVIEW OF THEORY 


The analytical development in the following sections is long and involved. Thus it 
seems worthwhile to present a simplified overview at this point with the aid of a series 
of figures. Figure 1 shows multiplication and time averaging of the governing equations 
for sound pressures p a and p^ to yield the p^pj^ correlation in schematic form 
(cf. refs. 3 and 4): the self -noise and shear-noise components of the source term 
(integrand) are singled out. Figure 2 brings out the principal new feature that arises 
when one microphone is replaced by two: the travel times from a source point to the two 
microphones are different in general. If the source region is noncompact (dimensions 
greater than a typical wavelength of the sound field), the difference in travel times will 
substantially reduce the p a Pj ) correlation. 

Figure 3 exhibits the phenomena in a more physical fashion. It shows how the 
source -observer difference in travel times r (which must be averaged over all source 
positions) gives rise to a cuspl in the correlation pattern as if/ — 0. The main point is 
that microphone b records somewhat the same pressure-time signature as microphone a, 
but with the signature delayed by an amount r that depends on the geometry. The cross- 
correlation p a p^ expressed as a function of time delay r is thus essentially the auto- 
correlation p a p a '( r ). Expressing t as a function of if/ ^and of 6 a , cp^, 0^ , cp^ as 
parameters) in effect maps P a P a ' ( T ) into p a p^ (if/, 0 a , etc.). In this mapping the cusp, 
which is characteristic of measured autocorrelations, reappears in the cross-correlation. 

The pressure signatures at microphones a and b due to a source element, although 
generally similar, differ by more than just the time shift r: they reflect the composite 
nature of the sources (fig. 4). The source strength density at a point consists of a super- 
position of six kinds of quadrupole T-j - three lateral and three longitudinal - some of 
whose directional patterns are shown in figure 4. At a given instant their relative ampli- 
tudes appear to be random; they are, however, related statistically, and the two- 
microphone correlation p a p^ involves a statistical average. 

The nature of the sources dictates the basic correlation that remains when the non- 
compactness effects are omitted: an example is the hourglass shape shown in figure 4. 

The basic correlation must be multiplied by the noncompactness factor to give the overall 
correlation. In the example case the forward lobe (if/ < 180°) will be large; the rearward 
lobe (if/ > 180°), very small. This will be alluded to in the Discussion, in the comparison 
of theory and experiment. 

*T he term cusp as used herein indicates a symmetric discontinuous finite slope 
of P a P b (^ at ip = 0, rather than a discontinuous infinite slope. 
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It may be of interest to return to the random (in time) distribution of six quadrupole 
directional patterns. In the present model these reduce to three ^ T ii> T i2’ T 13) * or the 
shear-noise contribution. Appendix A shows that T^g and T^g jointly produce a sound 
pattern of fixed shape (but fluctuating amplitude) whose directional lobes fluctuate stochas- 
tically in orientation about the jet axis. This is a variation of the "beam pattern" concept 
put forward by Maestrello and Pao (ref. 11), elaborated and refined in reference 12, and 
supported by experiment (ref. 13). The fluctuating beam pattern sweeps randomly over 
the two separated microphones a and b, and the time average of the contribution to p & p b 
is governed solely by the pattern shape and microphone positions. 

To round out the picture for the shear noise, refer now to the T 11 quadrupole: 
its beam pattern is of the form A(t)cos^0. This figure-eight beam pattern fluctuates 
solely in amplitude while maintaining a fixed direction alined with the jet axis. 

ANALYTICAL DEVELOPMENT 


Governing Equations 

Lighthill (refs. 8 and 9) has shown that the sound pressure radiated to a point x 
in the far field by a localized unsteady or turbulent flow is given by 


p(x,t) = (4 „c„2) 


tin 

at 2 


d 3 y 


(i) 


where T^ is a quadrupole strength density, 

T tj - "V] + T ij + ( P - c o 2 <>) 5 ij < 2 > 

that is normally dominated by the unsteady momentum flux P v i v y f° r example, in a tur- 
bulent jet at ambient temperature. Here t-j is the viscous compressive stress tensor; 
P, the local pressure; p, the density; c 0 , the ambient speed of sound; v^, the velocity; 
and 5.. = 0 or 1, as i ^ j or i = j. The brackets designate retarded time (to be spec- 
ified later); i,j = 1, 2, or 3; and repeated indices are summed over. The flow occupies 
only a limited region around the origin of coordinates. 

The quadrupole strength density T^ is approximated as P 0 v i v j and, following 
Proudman (ref. 10), x^v^x is observed to be the component v x of the resultant veloc- 
ity along x. For two points x a (at time t a ^ and x b (at time t^ - with 
I x a | = | x b | = x - this yields 
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( 3 ) 


P a ( = ( 4 to o 2x ) Po I 


Pb^ b ’^ = ( 47rc o 2x ) Po J 




8t 2 


k) : 


L at* 


,3— a 
d y 


d 3 y b 


where the brackets signify evaluation at the respective retarded times 


e a ,. -i*-n 

a a c Q 

f N b -y b l 
b % 


(4) 


Since x and x - y are nearly parallel for x in the acoustic far field, where x » y, 
these are approximately (fig. 5) 


t_ = 


a a c 0 c Q 


t - 1 - x i y b • 

b l b c Q + c Q 


(n a = x a /x) 
(n b = x b /x) 




(5) 


, o t lj 

Thus, the correlation of sound pressure at equidistant points x and x , at respective 
times t = t and t^ = t - r, is 


R(x a ,x b ,r) = p( x a ,t) p( x b ,t - t) 


2 , 


8t2 ( V *) . 8* 2 

o u + t 


,2 , v ,\2 


d 3 y a d 3 yb (6) 


—a -*-b 

The directions of the vectors x and x to the two observer points have been desig- 
nated by unit vectors n and n . The separation between the two source points y 
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and y b may be designated 77 and the midpoint by Y. The relations are 


— —a —b 

77 = y - y 


y = |(y a + y b ) 


( 7 ) 


—a — b 

The difference in travel times from the sources at y ,y to the respective observers 

T_ 

at x , x is, according to equations (5), 


—a -»a -»b -»b 

a = y • n _ y n 


_ Y • (n a - n b ) + 77 • (n a + n b ) 

2c„ 


( 8 ) 


If the time delay between receptions of the two signals at x and x is 


t & - t^ = r (specified) 


then the time delay between emissions from y and y is 


(9) 


*a “ t b (= r) = r + A 
Equation (6) may be transformed into 


( 10 ) 


R(jt a ,x b ,T) = p 0 2 (l6 A 0 4 x 2 ) A(v a ) 2 (v b ) 


d 3 77 d 3 Y 


( 11 ) 


2 2 A 

where is evaluated at y a ,t a and (v^j at y^t^. The flow is postul ated in the 

statistical sense as stationary and locally homogeneous; thus the correlation ^v^| ^v^j 
depends primarily on the time delay Tj and the separation 77. The brackets imply 
evaluation at a time delay = r + A. In this work the reference frame is stationary 
and the sources are stationary as well: they are specified as being unconvected by the 
mean flow. 

In what follows, a model for the two-point space-time velocity correlation 
( v x) 2 ( v x) 2 * s Postulated. The model is limited to the jet mixing region and varies 
with position Y in a fashion dictated by similarity concepts. A series of simplifying 
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assumptions, based on physical considerations, is applied to render the sixfold integra- 
tion tractable. With this same motivation the final analysis is restricted to zero time 
delay r between reception of signals at p a and p^. (This restriction may be relaxed 
without difficulty.) The final two-point sound pressure correlation thus reflects a number 
of approximations in the model. 

The instantaneous local velocity v is the resultant of the turbulence and stream 
velocities: 


v = u + U (12) 

with U taken parallel to the jet axis and u having zero mean. Since v x is the com- 
ponent of v along x, then 


v 2 = u 2 +2Uu + U 2 
x x x x x 


(13) 


(^) 2 K) 2 = K) 2 ( + + 


2 , \2 


. 2 / \2 


2 , „ .2 


■ 2 -a 


( u x) K) +K) «) +2u x( u x) u ; 


+ 2U x( U J) 2 “x + K <{<i + 2U x u x( u x) 2 


(14) 


Only the first two terms on the right-hand side of equation (14) will contribute to the jet 
noise. The three following terms are constant with r and their derivatives will vanish; 
the next two terms will likewise vanish, since the components of u have zero mean. 

The last two terms involve triple velocity correlations; it is argued in reference 4, 
appendix B, that their integral over 77 -space will vanish for homogeneous isotropic turbu- 
lence. These terms would also vanish prior to integration if the turbulence were postu- 
lated as a Gaussian process. 

On the assumption of normal joint probability for u^ and u^, equation (14) 
becomes (see, e.g., ref. 14) 
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I 


Self- Shear 

noise noise 



yfT T a TT b a b 
+ 4 U x U x Vx + 



+ Noncontributing terms 


( 15 ) 


The first term has no contribution from the mean flow U: this source term accounts 
for the "self -noise" of the turbulence alone. The second term ultimately depends on the 
mean flow shear in its contribution to the jet noise: it is termed the "shear noise." The 

third term is constant with r and will vanish on differentiation. 

Self -Noise 

The contribution of the self -noise to the correlation R(x a ,x b ,T) of sound pressure 
at the two equidistant points x 3 and x° with time delay t (the final evaluation spe- 
cializes to t = 0) is now formulated. The first term on the right-hand side of equa- 
tion (15) - the self -noise term - may be expanded as follows: 


a a a 
u x ~ n i u i> 


b a a 
u x = n j u. 


a b a b a b 

u u = n. n. u. u. 

XX 1 J 1 ] 


a b I b a b , a b a ,b 

u x u xj = 2 n i n j n k n Z u i u j u k u Z 


(16) 


However, the postulated isotropy of the u^ field may be exploited in order to obtain a 
simpler expression. To this end, new reference axes are chosen so that the 1-axis lies 
along x a , or n a , the orthogonal 2- and 3-axes being arbitrary: 







= i 



3 

= COS \p 

n 2 

- 


= 0 

b 

n 2 

- 

$ 

= sin \p cos p. \ 

n 3 


-a 

n 3 

= 0 

b 

n 3 


*b 

n 3 

= sin ip sin p. 



— — 




w — < 

J 


j 
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Sketch (a) 

Here (polar angle) is the angle between x a and x* 5 (or n a and n* 5 ) and p is 
an azimuth angle. (See sketch (a).) Then 


a b ~b a b ~b a b ~b a b 
u x u x = n l u l u l + n 2 u l u 2 + n 3 u l u 3 


(18) 


a b 
u u 
x x 


nl '-b'-b a b a b 

- 2 n i n j Vi u i u i , 


(19) 


The contribution of unit volume to the correlation (eq. (11)) to be evaluated with the 
use of equations (15) and (19) should be invariant with the angle p for fixed when 
convection effects are neglected: this is a consequence of the postulated isotropy of the 
turbulence. Thus there should be no change upon averaging over p. Performing this 
average on equation (19) with the use of equations (17) results in 


2 ( u x u x) =2 ( u l u l) C0S V + (v^) 


a b 
u l u 3 


sin^xp 


( 20 ) 


Following reference 4, it is now postulated that u^ Uj is factorable into a space 
factor appropriate to homogeneous isotropic turbulence (e.g., Batchelor, ref. 14) and a 
time factor: 


U i U j = S( T l) C ij W 


(21a) 



(21b) 
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(21c) 



f(? 7 ) = exp 



V 2 ='7 1 2 +1 2 2 +?7 3 2 


(2 Id) 


f'(v) = ~~ 

d? 7 


(21e) 


More generally, f may be any suitable universal function of ?? ; in the present model of 
turbulence it is taken to have the form shown, which has been used by Lilley (ref. 15). 

The longitudinal macroscale is denoted by L. 

It is known that the function exp^-o^iTjj j - which resembles the results of hot- 
wire measurements of turbulence - is unacceptable in the vicinity of Tj = 0; a zero 
slope is required there on theoretical grounds. However, experimentally the cusplike 
shape is observed extremely close to = 0; the region of failure is normally too small 
to be resolved from the measurements. Thus the specification (eqs. (21)) may be 
regarded as being applicable very, very close to, but not precisely at, = 0. The prac- 
tical effect of this exception for numerical work is negligible. (It may be argued that the 
foregoing implies a nearly singular behavior of 9%/a as — 0. However, restora- 
tion of the neglected time delay across a correlation volume should largely suppress the 
effects: t-^ can then vanish only at a single point within the correlation volume.) 

The self-noise contribution to the two-microphone correlation is obtained by 
inserting the first term of equation (15) into equation (11). Then with the use of equa- 
tions (20) and (21), there results 



x a ,x b ,r) = P 0 2 (A 0 V 



4 -2w f |r+A| 
e 


x 




( 22 ) 


where the retarded time difference r + A has been inserted for r^. 

An approximation may be made to decouple and simplify the rj and Y integrals. 
In the expression for A (eq. (8)), involving both rj and Y, the term in 77 accounts 
for phase incoherence in the source region. Suppose the coherent regions, or "eddies," 
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are taken to be acoustically compact (small compared with a typical wavelength of the 
radiated sound), then this term - in the postulated absence of convection - may be neg- 
lected. (An equivalent assumption has been made in the work of Lighthill (ref. 9) and 
others.) Equation (8) then reduces to 


A 



(23) 


In this approximation A is the difference in travel times of sound from a source point 
at Y — the centroid of an "eddy,'' or correlation volume - to field points at x a and 

i. 

x , respectively. The assumption of acoustic compactness, it is noted, is limited to the 
eddies: the total volume occupied by these eddies is not taken to be compact . 

Evaluation of the inner integrals in equation (22) covering all rj - space is elemen- 
tary, but a little tedious, and yields 


R S e(x a ,xV) = 


)3/2 7J .2c^4x2 


I + I COs2 * 




w 


t 3 4 

L u>^e 


• 2 "fl" a 1 d 3 ? 


(24) 


The remaining integral may be limited to the volume occupied by the jet, since u^ — 0 
outside the nominal jet boundary. 

It often suffices to avoid the integration in Y and look merely at the integrand in 
equation (24): this is the contribution to from unit volume of the jet. It was 

earlier remarked that the directional properties of the mean-square sound pressure 
emitted from a unit volume of a jet are fairly representative of those of the entire jet in 
a given (corresponding) narrow frequency band. However, no such simple correspondence 
is applicable to the two-point correlation of sound pressure, for reasons brought out 
below. Thus, the threefold integrals that largely sufficed for the single -microphone case 
must give way to a sixfold integral for the two- microphone case. Accordingly, the Y 
integration in equation (24) will have to be carried out. 

Suppose, just for the moment, that the major noise -producing region of the jet is 
postulated to be acoustically compact (and to contain the origin of coordinates) - that is, 
its dimensions are much less than the typical wavelength of the generated sound - then 
differences in travel times A are negligible in equation (24). In this case, the integral 
will reduce to a simple nondirectional factor. The factor (1/8 + 7/8 cos^if/) is thus the 
directivity of It for an acoustically compact jet. The integral, when A is not neg- 

lected, is therefore a multiplicative factor that allows for noncompactness of the jet. Note 
that it is not unreasonable for certain purposes to assume - as has been done here - 
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that the correlation volumes are compact; but it is totally unrealistic to assume that the 
many times larger jet is compact. 

Integral Over Noncompact Jet 
From equation (23) it can be seen that 


c q A a Y • n a - Y • n b 

measures the difference in distance from Y (center of an "eddy") to points x a and x* 3 . 
hence A measures the difference in travel times for sound originating from the eddy. 

This time difference reduces the cross -correlation of signals received at x a and 
x D in precisely the same fashion as time delay reduces the autocorrelation of a single 
signal (at x a , say) . This is apparent from the way A appears along with r in the 
exponential in equation (24). Thus the noncompactness factor can have a powerful effect 
in modifying the directional two-point correlation pattern predicted for a compact jet 
when A is large. 

Assume that the major contribution to R ge comes from the jet mixing region and 
limit the integral in equation (24) to that region. On abbreviating the integrand as g(y) 
and going over to cylindrical coordinates, the integral takes the form 



Gr dr dy dYj 


which can be transformed into 


(25) 



GrYj dY x dy d 



(26) 


The scale L and characteristic radian frequency depend essentially on Yj alone, 
and 


Ui 2 (y) = Ui 2 (R) F^j (27) 

describes the turbulent velocity in the spreading annular mixing region: it specifies sim- 

ilar radial profiles at successive axial distances. To the present accuracy the part of the 
integral in equation (24) varying with r is just 
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( 28 ) 



Since F peaks strongly at r = R, this is approximately Rgp where g^ is a pure 
number. Additionally, although and L vary with Y^, their product is essentially 
invariant (constant Strouhal number). 

Thus the correlation (eq. (24)) reduces to 


R 




) R (uj 2 (R)] /j [ ? 


2 4 2 

77 C Q X 


(| + | c °sV) 


— dY ! f 

Jq L 1 Jf 


277 -2u>j|t+A| 


dy 


(29) 


where £ is the effective length of the source region. An approximation to the mean of 
diverse hot-wire measurements is the linear dependence L = Y^y^, where gg is a 
constant of order 5 to 10. Accordingly, the correlation may be written as 



-a — b \ 
x ,x ,t) = 



(I + 2 cosV) 


xif'dY 1 f 2 ' e - 2 “f M 
£ -'o 1 277 Jo 


dy 


(30) 


For evaluation, a suitable expression for the difference in travel times A is 

—►cl t To 

required. Spherical polar coordinates are introduced for n and n and the cylindri- 
cal coordinates for Y are recalled: 


Y 1 = Y 1 


Y 9 = r cos y 


Y„ = r sin y 


x a ^ 

a X 1 

m = — = cos 0 Q 
lx a 


x a 

n 2 = lT =sin 0 a cos ^a 


x a 
a X S 

n„ = — = sin sin c>„ 
o x a ^a j 


(31) 


with similar expressions for ^n^, ng, n^. 
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It has already been noted that F(r) in equation (27) peaks very sharply at r = R. 
This implies that the sources of sound are concentrated in the neighborhood of the cylinder 
r = R. The expression for A is therefore simplified with very little error in final 
results by replacing r by R in equations (31). Then A may be expanded as 

c q A = Y • ( n a - n b ) = Y-^cos 0 & - cos 0 b ) 


+ R cos 
+ R sin 


y^sin 0 a cos <p a - sin 0 fe 
y^sin 0 a sin <p a - sin 0 fe 


cos <p b ) 



(32) 


This reduces to different forms for the two special cases under consideration (fig. 6): 
(1) the two microphones lie on a circle of latitude 


lie on a meridian circle 
the jet axis. 


0 a = 0^ and (2) the two microphones 


cp a = (p b - Note that the polar angles 


Vt 


are measured from 


The double integral in equation (30) is evaluated in appendixes B and C for cases I 
and n, respectively, by using certain simplifying assumptions. The evaluation is limited 
to the special case r = 0. (That is, the time delay of signal reception between the two 
microphones is specified to be zero.) Insertion of the final results into equation (30) 
gives 


Case I: 0 & = 0 b ; <p a = 0 



a — b 


o) 


= g l g 2 P Q 2 / CtJ f L '\ 4 
2^/2 7T \ C 0 / 



(W c°sV)l(q S e) 


•\ 




sin e 

de 


Noncompactness factor 




(33) 


%e - 4 <"fX B / C °) 


sin 0. sin 
b 


^bl 
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Case BE: <p a = cp^ = 0 

»(*.*>. o) - 

sc' ’ > ' 2 l/2 it ^c„ 






R£ 


x 


[H c ° s2 ( 9 b - «»j) (l ~ 2 e 6 ' 26) 


Noncompactness 

factor 


where 


6 = 


( w f> 


COS 0^ - COS © a 


(34) 


In both cases signifies an effective average of coj over the length It of 

the mixing region (appendix B). Equations (34) become invalid near 0^ ~ -0 & and for 
that subcase are to be replaced by equations (33) with 0^ = 9 , cp-^ = 180°. 


Shear Noise 

The shear-noise contribution to the two -microphone correlation is obtained by 


replacing ^v a 

tion (15), namely, 4U a U b u a u b Again, reference 4 is followed in describing the mean 

X X X X 

flow shear by the correlation 


J / b\ 2 

(x ) * n ec l ua ti on (11) by the second term of its expansion in equa- 


U 





(35) 


where the prime designates the second point, rj n is measured radially outward, and the 
flow is treated as locally planar within a correlation volume; correspondingly, 


u a u b 

x x 



(36) 


n b b cL b sl b t \ 

For u u = n. n. u, u, , it is again specified that u. m has the form C..(tj) g(T), which 

XXlJlj 1J lj'' 

is spelled out more fully in equations (21). Putting this all together gives the source term 
as 
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(37) 


^3 " 


, ab ab 
4n^n^n^ nj 




This is to be inserted in equation (11) and evaluated at Tj = r + A to yield the two- 
microphone correlation. Specializing to the simplest case, r = 0, gives 



ab„a„b n 2 

n i n i n i n j pQ 

4ff2c 0 4 x2 



x 



d 3 77 d 3 Y 


(38) 


For reasons related to the axisymmetry, cylindrical coordinates If = Y.., r, y will 

o— ► 

now be used to locate the volume element d°Y. The integral in 77 -space will be simpli- 
fied if the coordinate axes are rotated counterclockwise through angle y [77" ~j)\ 
then the radial two-point separation rj n will become but the 0^(77] will be 
unchanged in form because of the postulated isotropy of the turbulence. All but two of 
the terms will have zero integrals; the surviving integrals are 



L 3 o 


2(1 + a) 3 / 2 


(39) 


Upon insertion of equations (31) referred to the rotated axes {cp — along with equa- 
tion (39), and noting that <*&()[) vanishes outside the jet, equation (38) reduces to 



a 


877^C 0 4 X^(1 


a)3/2 Jjet 


c ( 


2 2 
cos 6 a cos 0^ 


+ i cos 
2 


0„ cos 
a 


sin 0 o sin 
a 


x sin (p sin 


- \ 2 T 3 4 

cpb) L co f 


-w f |A| 


U‘ 


: (y) 


d 3 Y 


(40) 


Reversion to the original reference axes gives, on noting that <p and cp^ are angular 
deviations from the angle y orienting the volume element, 
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R sh (x a ,x b ,°) = 


87T 2 C 4 X 2 (1 


. + a) 3/2 ^jet C 


2 2 1 

cos G cos + - cos G n cos 0, sin A sin 0, 
a b2 a b a b 


X sin(<? a - y) sin^ b - y) 


2 T 3.. 4 " w f| A | tt 2/- 


Uj L aj f e 


U 2 (y) d 3 Y 


( 41 ) 


Just as for the self -noise let 


i* t d3? - OT Jet GrY i dy i ^ d (iff 


(eq. (26)) and invoke the same similarity assumptions concerning L, u^, and u^ 
leading to equation (30). Additionally, the specification <P a = 0 can be made without 
loss of generality: <p b will take on the role of <p b - cp a because of the axisymmetry. 

With these stipulations together with u^ 2 = u^ 2 (R) F^r-R^Y^j and U 2 = Constant 
(annular mixing region), equation (41) reduces after the radial integration to 


R 


sh' 


(■* a,b„) M 4 r< uS ! w. lf i r 2l, C 

'' ’ ’ 8» 2 c_ 4 x 2 U + <t)3/2 l j 0 ‘Jo L 


cos 2 0„ cos 2 0, 
a b 


+ cos # a cos 0 b sin sin 0 b sin 


(y - cp b ) sin ^ 


■co f |A 


d y 


(42) 


where 


gg = Yy/iL = Number (scaling law) 

go = T . f(^\ = Number (since F = 0 at lower limit) 

3 J-r/Yj \ y i/ \ y i) 

and A^y; is defined in equation (32). 

The experimental evidence indicates that the compactness condition 
^u>j|A|s47fR/X«l^ is not normally met. Thus the exponential in equation (42) that 
accounts for noncompactness of the source region (cf. earlier discussion regarding self- 
noise) must be retained. The integration is discussed in appendix D. The rather lengthy 
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results for the shear -noise correlation R g k(x a , x b ,o) will not be repeated here. 
Instead, the combined contributions of self -noise plus shear noise will be summarized in 
the following section. 


Summarized Predictions 


In order to put the final equations in the most compact form, a normalization pro- 
cedure is used. It is observed that the correlation R(x a , x b ) = PgP^ becomes merely 
p 2 (x a ) when x b = x a . On specializing to the self-noise, the correlation in equations (33) 
reduces to 



(43) 


-♦■fa 1 3, 

A similar reduction occurs for the shear -noise correlation when x = x ; the maximum 
value occurs with x a oriented along the jet axis (0 = 0): 



a —a 


0 ) 


max 



max 


'bijh t L \ 4 u i 2 u 2 rip 

ilt \ c o / X^(l + cr)^/^ 


(44) 


The ratio of the peak shear-noise intensity to the self-noise intensity is a pure number; 


£ = ( Psh lnax = 2 -3/2 ^3 a 

2 d 2 g l 2 (1 + a) 3 / 2 

p se 1 


(45) 


When the integrals gj and gg over the jet radius are replaced by unity, equation (45) 
reduces to the corresponding ratio B/A referred to unit volume given in reference 4. 
The parameter cr in equation (45) is evaluated in terms of the mean flow shear 3U/9r 
in appendix B of reference 5. 

The final equations for the combined contributions of self -noise and shear noise to 
the two -microphone correlation 


R(x a ,x b ,o) = R ge (x a , x b ,0) + R gh (x a , x b ,o) 


are expressed in terms of pjL. 
are 


For the two cases of main interest (fig. 6), the equations 
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Case I: 9 & = 6> b = 9, cp^ = 0 


(x_^c ? o) = (l + 7 cos 2 V/)l(q se ) + |(cos 4 0 + cos 2 0 cos ^)l(q gh ) 
P se 

9 sin 2 0[l(q sh ) - J(q gh )j 




- ^ cos 2 ' 


where 


I(q) = - r /2 e - q Sine de 
71 Jq 

J(q) = i f’ /2 U s ),-<1 sl " « d < 

7T Jq ^ ' 


■\ 


/ Ncncompactness factors 


> (46) 


q sel 4 


q c u i 2 


l <c0f>(R/c o ) 


l sh 


. ^ b 

sin 6 sin — 
z 


= Effective average of over mixing region 


-*>a -►b 

= Angle between x and x 


= cos 


- ^cos^0 + sin^0 cos <p^j 


J 


The noncompactness integrals I(q) and J(q) are functions of 6 and cp^ that are 

readily evaluated numerically by computer. For fixed 0 they approximate roughly to 
-const \cp I 

the form e I b | ? the constant depending on 6 . 

Case H: <p h = cp a 

= 1-^(1 ♦ 1 c 0S 2 10 b - O 
T" 26 8 8 \ b a) 

n — 


R 


se 


B 1 - e“^ / 2 2 1 

+ ^ — (cos 0 a cos ^ + 2 cos cos ^b s * n ^a s * n ® 


*>) 


where 


(47a) 
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(47b) 


6 = 



- cos 6 

ct 


(Eqs. (47) are invalid in the vicinity of 
with = 180°; cf. footnote 2, p. 37.) 

(l - e" 26 )/26 and (l - e" 5 )/6; they 
x b — x a ^or as 0^ — 0 . 


0^ = - 8 & and are to be replaced there by eqs. (46) 
The noncompactness factors in this case are 

again have an angular dependence with a cusp as 


DISCUSSION 

Comparison With Experiments of Maestrello 
The theoretical equations for the two-microphone correlation 

P a (x a ,t)pjj(x b ,t) 55 R(x a ,x b ,C>) 

are given in equations (46) and (47), with ancillary definitions in equations (43) to (45). 

The correlations are nondimensionalized by division by the mean-square self-noise 
o 

sound pressure p to simplify the expressions. However, it is preferable to use the 
more usual scheme of normalization: 



(48) 


The values of R may be taken from equations (46) and (47) as the case may be, and it 
2 

will be noted that p will cancel out. 

SG 


The normalization (eq. (48)) corrects for large amplitude changes resulting from 
convection and refraction effects excluded from the theory. To the extent that they could 
be approximated as uncorrelated multiplicative factors for p a and p^ herein, they 

would cancel out along with p^ . 

SC 


The governing parameters in the equations are 2 for case I = 0 & = 0) 

and for case H = <P a y In terms of Strouhal number Ng tr sf^iyUj, these 

take the respective forms 27 rMj^Ng tr ^ for case I and 27rM^(f/D)^Ng^ r ^ for case n. 

Here D is the jet diameter, £ is the effective length of the noise -generating "cylinder," 
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and the angular brackets indicate an effective average over length £. The ratio f/D 
was taken as 5.25, about the length of the mixing region. Mach number enters indirectly 
through its influence on {Ng^ as well as directly; the effects may be opposing. For 
numerical application of the theory at Mj = 0.75, was taken as 0.40. This was 

frankly chosen for best fit to the experimental data, but the value is compatible with tur- 
bulence correlation data. 

The turbulence data of reference 16 may be drawn on to elaborate the last remark. 
An appropriate autocorrelation is one referred to a reference frame moving with the local 
convection velocity. This moving -frame autocorrelation is the envelope of two-point 
correlations measured in the fixed (laboratory) frame (ref. 16). These envelope auto- 
correlations in reference 16 and elsewhere roughly resemble the simple exponential 
choice for g(T) herein, but they deviate at small r and show a sharp upturn as r — 0. 
Thus the experimental correlations are better fitted by the sum of two exponentials: one 
with fast decay to dominate the slope and higher derivatives for small r plus another to 
dominate at large r. 

Now the key features of the theory are displayed at small r (including the cusp 
at t = 0); hence, since the theory (in its present version) is limited to a single exponen- 
tial, it is appropriate to choose the exponent for best fit to the experimental moving- 
frame autocorrelations in the region of small r. This has been done with the envelope 

curve of figure 21a of reference 16, which is a plot of UjU^Uj^/D^ for rj = 0 in the 

present notation; the measurements are made at = 4.5 and r/R = 1.0, with 

M.j ~ 0.13. The fit is made to the dashed portion of the curve at small r terminating 

o V\ “CfJr j 7" I 

in the point Lby D = 0.065, u \ u i = 0-871; thus e 1 is chosen to pass through this 
terminal point. The result, with fj = is 

fjD/Uj - N slr = 0.337 
This is to be compared with 

= {^str^ = ®*^00 ( Y l/° ay e ra ged from 0 to 5.25^ (50) 

used in the theory. Since Ng^ r increases with decreasing Y^^D, the two values are 
quite compatible. 

Values of f^Dy^Uj quoted in the literature (cf., e.g., ref. 3, p. 120) are obtained in 
a fashion that is different and not appropriate to the present choice for the theory. They 


(Yj/D-4.5) (49) 
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do not result from a fit to the experimental autocorrelation at small r. Instead, 

= 27 rfj is chosen as the reciprocal of the time -scale: the integral of the autocorrela- 
tion. This yields manyfold -s/raZ/er values of f*. 

Experimental measurements of the two-microphone correlations are available from 
a massive and definitive investigation by Maestrello (ref. 7). Only a few of his results 
are drawn on for the comparison. In particular, attention is limited to the peak values of 
the correlations when plotted against time delay t. In the predictions of the theory 
these peaks occur at r = 0. In the experimental measurements the peaks are essentially 
at r = 0 for the axisymmetric case I (explained in fig. 6), but are somewhat offset from 
r = 0 for the nonsymmetric case II (fig. 6), owing to convection and refraction effects 
excluded from the theory. 


Refer once again to figure 6 to see clearly the two main cases of microphone 
positioning. In case I both microphone position vectors x a (to p a ) and X* 3 (to p^j 
lie on a cone (0 & = 0^ = 0; <p a = 0, <p b varies). In case n, the vectors x a (to p & ^ 
and x b (to p^) lie in the plane of the jet and p a and p^ move around a meridian 


circle (p a at 0 & , ^ at 0 fe ). 


Example polar plots of Maestrello' s experimental measurements (ref. 7) are shown 
in figure 7 for case I (0 & = 6^ = 0; <p a = 0, <p^ varies) at Mach number Mj = 0.75: the 

normalized two -microphone correlation is plotted against the azimuthal separation cp^ 
of the microphones. These and other cases are replotted in Cartesian format and com- 
pared with predictions of the present theory in figures 8 to 11. The general shape of the 
experimental curves, in particular the cusp at cp^ = 0, is well predicted. Quantitatively 
the agreement is reasonably good in the downstream quadrant (fig. 8), excellent in the 
90° plane (fig. 9), and only fair in the upstream quadrant (fig. 10). It is significant that 
the theory, which omits convection and refraction effects, is most accurate in the 
90° plane, where such effects should be minimal. 

An additional factor presumably contributing slightly to the accuracy at 90° is the 
vanishing there of the slightly uncertain amount of shear-noise contribution. In the speci- 
fications herein for the theory there is a parameter /3/2 = Peak shear-noise intensity/ 
Self-noise intensity. Values of /3 inferred from single-microphone measurements of 
jet noise vary from 1.58 to 3.55 (ref. 5); the theoretically based value /3 = 2.0 was used 
herein, except where otherwise specified. 


To elucidate the significance of the shear-noise/self -noise ratio /S/2 further, the 
effects of three different values were explored: theoretical curves for /3 = 0, 2.0, 

and 3.55 were evaluated for case I with 0 = 30° and 60°. These are presented in fig-' 
ure 11, along with Maestrello's measurements. The choices /3 = 2.0 and /3 = 3.55 
give comparably good agreement with experimental results at both 0 = 30° and 0 = 60°. 
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On the other hand, the choice /3 = 0 (no shear noise) considerably degrades the agree- 
ment at 9 = 30°. Thus the shear noise cannot be ignored, but the results are not highly 
sensitive to its precise proportion, except that values of /3 much below 2 must be ruled 
out. Put another way, the results are not highly sensitive to the basic directivity 
resulting from the "nature of the sources" referred to in figure 4. 


The agreement between the theoretical curves and the experimental values in the 
upstream quadrant 9 > 90° (fig. 10) is substantially poorer than in the downstream 
quadrant 9 < 90° (fig. 8). In particular, the progression with increasing 9 is wrong. 
This can be explained qualitatively in terms of convection and flow- shrouding effects 
omitted from the theory. A crude accounting for these Doppler effects predicts changes 
in the effective difference in travel times A to the two microphones: it seems that a 
corrected A would be smaller than that of equations (35£) for 9 < 90° and larger for 
9 > 90°. Such changes are in the direction to improve agreement with experiment. 

Polar plots of Maestrello's experimental correlations for case n = <p a ; 
f? a fixed, 0^ varies) are shown in figure 12 ^Mj = 0.75). These are replotted in 
Cartesian format and compared with theory in figure 13. The agreement for 9„ = 60° 

a 

and 120° is seen to be excellent for moderate values of |0 U - 6 I. For 0 = 30° the 

o o I b a l a 

agreement is again excellent in the range 30 s < 60 , but the theory breaks down 
badly in the range -30° = < 30°. The breakdown here is not unexpected, since this 

range defines the M refraction valley"; within this valley refraction effects, excluded from 
the theory, become dominant as 0^ -*■ 0. 

More explicitly, within the range -20° < 6^ < 20° the experimental correlations 
show a deep dip unmatched by the theoretical curves. The loss in correlation is thus 
associated with those quasi-axial sound rays that have long paths through the jet turbu- 
lence. It may be speculated that phase modulation by the turbulence (that is, line 
broadening) accounts for the extra loss of coherence represented by this dip. 


Maestrello's measured correlations at Mach numbers of 0.59, 0.88, and 1.00 
(ref. 7) were broadly similar to those at Mj = 0.75, but with substantial individual varia- 
tions. No quantitative attempt was made to apply the theory at these other Mach numbers. 
Since Mj is not a direct parameter of the theory, the only significant free parameter 
that may reflect an influence of Mj is the Strouhal number Thus one could 

vary ^f^Dy^Uj^ from the value 0.40 used for Mj = 0.75 to different values in attempting 

to make calculations for different Mach numbers. 


Qualitatively it would seem that the agreement with Maestrello f s results, with 
<f f D/Uj> chosen for best fit, will be degraded at the higher Mach numbers. It is the 


author’s view that convection, refraction, and diffraction effects, omitted from the theory, 
are largely responsible for this; these effects, of course, increase progressively with 
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Mach number. Another view has been suggested by Maestrello. He attributes the devia- 
tions with Mach number to the term in rj in the time-delay equation (eq. (8)); that term - 
neglected herein for mathematical simplicity - accounts for sound travel time across an 
"eddy." 


Further Considerations 

The theory developed herein is not exact: it embodies a number of simplifying 
assumptions. The extent of agreement with experimental results must therefore be 
assessed in the light of these approximations. It is worth taking a closer look at some of 
the approximations at this point; to this end, the statement of assumptions in the 
Introduction is amplified in the following paragraphs. 

Flow refractive effects are suppressed because the governing equation - the 
Lighthill equation - is not a convected-wave equation. Convection of the sources is 
likewise not accounted for. The major contribution of these effects is thought to cancel 

on nondimensionalizing by \j p & 2 p^. That is, although convection and other flow 

effects can powerfully affect p a and p^, these effects are expected to have rather sim- 
ilar multiplicative effects in the numerator and denominator. 

The multipliers resulting from flow effects will not, however, cancel exactly. Put 
another way, the flow effects will alter the predicted non dimensional correlation some- 
what. This comes about because both the source and sound wave convection by the flow 
will alter the effective difference in travel times A to the two microphones. As men- 
tioned earlier, a rough assessment of these changes seems to be in the direction to 
improve the agreement with experiment. 

Another principal limitation is the modeling of the two-point space -time correla- 
3 . b 

tions of turbulence u j u j • These govern the character of the source terms in the 

Lighthill equation. The functional form assumed in equations (21) embodies certain sim- 
plifying assumptions: chiefly the assumption of isotropic turbulence and a specific fashion 

g( t) in which time delay r enters. The selection of g(r) as a multiplicative factor 
greatly simplifies the mathematics, but it is bound to be an oversimplification. And the 
choice of g(r) as a simple exponential decay in equations (21) approximates the general 

cl b 

form of the hot-wire measurements of u. uV (t) autocorrelations (referred to the con- 
vected reference frame), but not the detailed shape, as discussed earlier. 

-CO* I T I 

The defects in the chosen g(T), which is e 1 , become apparent when the 

microphones at x and x are brought together; then R(x ,x ,t) becomes the auto- 
correlation R (x a ,x a ,r). The time difference A in equation (30) vanishes and 
-2 { au? \ I r] 

e v l/ - an average over the jet - comes outside the integral; this will define the 
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shape of the self -noise autocorrelation. Similarly, one would obtain e 11 for the 
shear -noise autocorrelation. (In the general case, the predicted autocorrelation is pro- 
portional to d^g^(T)/dr^ for the self -noise and to 8 %(t)/3t^ for the shear noise.) 

A typical measured autocorrelation from Maestrello's investigation (ref. 7) is the 

curve labeled cp^ = 0° in figure 14. It does resemble e ( co f ^ I T I for moderate t, 

but for large t it has a negative loop. Working backward from this, it could be argued 

that g(r) should be chosen so that the derived autocorrelation ^ which for 9 ^ 90° is a 

composite of d^g jb-A and b^g^j bT^j should match the experimental autocorrelation in 

having a negative loop. The new g(T) then becomes g(r + A) in the formulas of the 

offer' 

theory. With this replacement the negative excursions of figures 7 fPjurts would seem to 
have an explanation: the negative loops in g(-r) would map into the negative values of 
P a P^ as A maps into <p^. 

The scenario of the last paragraph is, however, oversimplified. The time factor 
g(r) refers to unit volume of the jet turbulence, whereas the measured autocorrelation 
refers to an average over the jet. The spectrum from unit volume - of which b^g(r)/bT^ 
is the Fourier transform - is narrow band, whereas that from the entire jet is a broad- 
band summation of these narrow bands. The autocorrelation, being a Fourier transform 
of this broad-band spectrum, will differ in character from d^g/br^. Thus the corre- 
spondence relating an average b^g/br^ to the experimental correlation alluded to in the 
last paragraph is not justified in a detailed way. Nevertheless, the final sentence of the 
last paragraph is still believed to provide an explanation of the negative excursions in 
figures 8 and 11, but in a general, rather than detailed, fashion. 

Characterization of the source region as noncompact and essentially incoherent 
neither proves nor disproves a significant role for large-scale coherent structures in the 
generation process. Spatially, the reported structures are more long-scale than large- 
scale: they are characterized by long, convoluted vortex filaments. When included in 
more random turbulence, the structures do not appreciably modify measurements of 
instantaneous two-point correlations; such measurements in a jet indicate spatially 
coherent regions of only modest dimensions, for example, less than the jet diameter in 
the mixing region. Thus the effective average size of the coherent regions is small for 
sound -gene ration purposes. 

Another interpretation is this. A large-scale coherent structure may display a 
certain time scale when sensed by pairs of probes of fixed separation following the mean 
flow, but when sensed at points each of which follows the local flow, will exhibit a stretched- 
out or extended time scale. It is in this time sense that the structures are especially 
coherent. Now although the large spatial scale (say L') of the coherent region increases 
its radiative efficiency, this will be more than offset by its extended time scale or reduced 
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characteristic frequency (say co f '). In fact, the mean-square pressure at afield point 
emitted from unit volume is proportional to (e.g., ref. 17). Thus the large- 

scale coherent structures are not particularly effective sound generators. An exception 
to this statement might be made for axisymmetric structures, for example, ring vortices; 
however, Maestrello' s measurements are incompatible with a ny a ppreciable component 
of axisymmetric emissio n in the jet noise . 

CONCLUSIONS 

An extensive series of measurements of two-point far -field correlations of jet 
noise have been carried out by Maestrello and reported in NASA TM X-72835. The 
results displayed some surprising features, and this motivated the present attempt at 
a theoretical explanation. The predictions of the theory have been compared with 
Maestrello' s measurements, and the following conclusions are drawn: 

1. Despite a number of simplifying assumptions the predictions of the theoretical 

model conform well with experimental results: (a) they exhibit the characteristic cusp 

in the curves of two-point correlation plotted against microphone separation, and 

(b) more than that, in most cases the quantitative curves are predicted reasonably well. 

2. The analysis indicates that the shape of the curves is dominated by differences in 
travel times from a source to the two microphones. This is a consequence of noncom- 
pactness of the source region compared with a typical wavelength of sound. It shows up 
most markedly in the relative narrowness of directional lobes in the plane of the jet axis 
(case II) as compared with those in planes perpendicular to the jet axis (case I). These 
features are correctly predicted by the theory. 

3. Characterization of the source region as noncompact and essentially incoherent 
neither proves nor disproves a role for large-scale coherent structures in the noise- 
generation process: the reasons are detailed in the Discussion. However, the demon- 
strated absence of axisymmetry in the correlation pattern rules out any significant role 
for axisymmetric source structures. 

4. As a further corollary to conclusion 2, the basic directivity arising from the 
composite nature of the sources makes but a minor contribution to the correlation curves. 
It is, however, distinctly discernible in the curves of case I at a polar angle of 90°. 

5. A variation of the "beam pattern" concept of Maestrello and Pao (J. Acoust. Soc. 
America, vol. 57, no. 4, Apr. 1975, pp. 959-960; NASA TN D-8104; and NASA TN D-7269) 
accounts for the basic directivity. The sound emission from each small volume is seen 

to consist of a rather directional quadrupole pattern whose lobes (‘ ' fluctuate 

ar\d 3 (rip ///> u.ac. 

randomly in direction. Computation of the time average recovers the smooth patterns of 
basic directivity obtained herein by other methods. 
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6. Effects of convection of the sources and the sound waves by the flow are excluded 
from the theory; they are expected to cancel partially in the normalization of the correla- 
tions. There will be a residual effect resulting from altered travel time differences to 
the two microphones. This is estimated to be in the direction to improve the agreement 
with experiment; however, its inclusion would pose a formidable task. 

7. The foregoing conclusions refer to the two-point correlation at zero time delay . 

An alternative format - with time delay reintroduced - is its Fourier transform, the 
cross -s pectral density ; this separates out the contributions of different frequencies. The 
differential travel time effect (cf. conclusion 2) would be expected to be relatively unim- 
portant at low frequencies (e.g., at wavelengths equal to or greater than 2ti times the jet 
diameter). Maestrello has recently made Fourier transforms (as yet unpublished) of his 
correlation measurements which indeed appear to show this. Thus, conclusions 2 and 4 
do not carry over to the cross-spectral densities at the lower frequencies. 

By way of closure a general observation may be made. A valid physical model of 
jet-noise generation must be compatible with the major features of both single -microphone 
intensity patterns and two -microphone correlation patterns. The model central to the 
present work, developed in the sixties, meets this criterion. The newer models in the 
literature remain to be tested. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
September 13, 1976 
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APPENDIX A 


QUADRUPOLES WITH FLUCTUATING ORIENTATION: 

THE BEAM PATTERN CONCEPT 

The beam pattern concept discussed in the body of the paper is most simply 
illustrated via the shear-noise quadrupoles. According to reference 4, equation (15), 
the contributing quadrupoles are 

T 11 = 1 T 12 = P v l v 2 T 13 = P v l v 3 

= pUUj = pUu 2 = pUu 3 (Al) 

For a general quadrupole T.j occupying unit volume, the contribution to far -field pres- 
sure is 

x.x. 9 2 T,, 9 2 T.. 

p 5i rc — — i ^ - Directional factor x JJ- (A2) 

i x^ 9t^ 9t 2 


where 


Thus 


x^ = x cos 9 

Xg = x sin 9 cos cp S 

Xg = x sin 9 sin cp 


x l x l 

x l x 2 

x 2 

x l x 3 

x 2 


cos 9 


cos 6 sin Q cos cp > 


cos 9 sin 9 sin cp 

J 


(A3) 


(A4) 
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where 


p gb a aj(t)cos^0 + agWcos 9 sin 9 cos q> + a 3 (t)cos 9 sin 9 sin <p 
a x (t) = T 11 

& 2 (t) = t 12 
a 3 (t) = t 13 

" indicates double time differentiation. 

This may be rewritten as 

p gb cc a^(t)cos^0 + Aj(t)cos 9 sin 9 cos (<p - e(t)) "'j 




tan e(t) = 


a 3 (t) 

a 2 (t) 





(A 5) 


(A6) 


The first term accounts for the T^ quadrupole: its cos 2 0 pattern merely fluctuates 

in amplitude. The second term accounts for the Tj 3 and T^ 3 quadrupoles jointly: it 
has the directional pattern of a T 12 quadrupole in a reference frame rotated by a fluc- 
tuating angle e(t). In effect the beam pattern of this Tj 2 quadrupole is fluctuating in 
direction somewhat as pointed out by Pao and Maestrello (refs. 11 and 12). 

Sketches of the directional polars of the two components of the instantaneous shear- 
noise pressure p gb are given in figure 15. 

The two-point correlation of far -field pressure is now formed according to this 
fluctuating beam pattern model. From equations (A5) with x a °c <p a ; x b cc 6^,^; 
for case I (0 a = e b = e, <p a = 0): 

Psh^^Psh^ 1 ^ a j 2 (t)cos 4 0 + A 1 2 (t)cos 2 0 sin 2 0 cos (<p b - e(t)]cos e(t) 

+ Cross-product term 

cc a 1 2 (t)cos 4 0 + | A 1 2 (t)cos 2 6> sin 2 e|cos ^ + cos - 2e(t) 

+ Cross-product term (A7) 

The average P sh (x a )p sh ( ^ B R s h(x a ,x b ,°) is required. The isotropy of the 
turbulence dictates that all angles e are equally probable irrespective of the 
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o 

amplitude Aj . Thus, the average over the amplitudes may be carried out first - so 
2 — *7 2 2 

that a^ (t) — a^ , A-^ (t) — Aj - before carrying out the average over e. The prob- 
ability that e lies in the range de is de/2;r, so that 


R gll ( x a , x b ,o) °c a^ 2 cos 4 0 + i Aj 2 cos 2 0 sin 2 0 cos cp ^ 

+ I Aj 2 cos 2 e siA C 2 ' cos («, b - 2e) || 


+ Integrals of cross-product terms 


(A8) 


2 2 

The integrals over e all vanish; also it is clear from reference 4 that a^ and A^ 
are equal. Hence, finally 



4 12 2 

cc COS 6 + - C0S 0 sin 0 COS 


(A9) 


Alternatively, the correlation can be expressed in terms of the angle \J/ between the 
—a — b 

vectors x and x , by use of equation (D5), in the form 


4 2 

td /V.a ~b ^ cos 0 + cos Q cos 
R sh' x > x 00 5“ 


(A10) 


Equations (A9) and (A10) re-cover results obtained in the body of the paper by a 
more conventional procedure. Thus interpretation of the source terms of jet noise as 
quadrupoles with fluctuating orientation - the fluctuating beam pattern concept - is 
given strong support. 
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APPENDIX B 


EVALUATION OF NONCOMPACTNESS FACTORS FOR CASE I: 

e a =e b = 0 > *a“° 

The double integral in equation (30) defines the n oncompactness factor for self- 
noise. The inner integral, to be evaluated for r = 0, is 


Self -noise 


1 P 277 -2o> f |Aj 
I„„ = — \ e d y 

2n Jq 


se 


Closely related integrals must be evaluated for the case of shear noise: 
Shear noise 

1 p 2ir — ol»£| A. | 


sh 


~ j_ r 

2ir Jq 


dy 


, ^27T , -ojJaI 

sln ( y ‘f) e dr 

For the present case 0 a = = 9, cp a = 0, equation (32) for A reduces to 

cos y (cos - 1) + sin y sin 


A = - — sin 6 
c o 


— sin 9 sinl 
c o 


ln (?) sin ( v * 


^b 

2 


Put 



sin 9 sin| 




e = y - 




(Bl) 


(B2) 

(B3) 


(B4) 


(B5) 


(B6) 
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Because of the tt/ 2 periodicity of the integrand, the limits of integration can be divided 
by 4 as well as arbitrarily shifted in phase. The phase shift can be chosen so that A is 
always positive and the | j can be dropped. The integrals can then, with the use of 
equations (B5) and (B6), be written in the more convenient forms 


o r n/2 -q sine 

- *0*— ) “ I Jo e de 

2 r n/2 -q .sine 


T Tin \ 4 f ^ 2 2 P / q sh Sin6 . 

J sh " J ( q sh) " i J n S1 " fe de 


(B7) 

(B8) 

(B9) 


Although these are very suggestive of defining integrals for Bessel functions, the limits 
are different and no such reduction appears possible. The integrals are, however, very 
amenable to numerical evaluation. 

Equation (30) requires that I (eq. (B7)) have a further integration over Y^; 
there are similar requirements for 1^ and J^. If (to which q gg and q g ^ 
are proportional) is taken as a constant independent of Y^, the integral merely yields a 
factor of unity. For a jet, however, it would be more realistic to take o>j cc Y^“ 4 
(cf. ff. eq. (28)), but the price would be a more troublesome Yj -integration. This extra 
complication was not considered warranted in view of the other approximations in the 
analysis. Instead, oj ^ is written as a constant in equation (B5), and (Wf) is 

called an "effective" Yj-average. 
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EVALUATION OF NONCOMPACTNESS FACTORS FOR CASE H: 

<?b = <?a 

The double integral in equation (30) defines the noncompactness factor for self- 
noise. The difference from the case treated in appendix B is in the angle specification 
governing A. There will be no loss in generality on replacing ^ = <p a by ^ = (p & = 0, 
because of the axisymmetry of the jet noise. Then equation (32) for A reduces to 


c q A = Yj(cos - cos + R cos y (sin 0 & - sin 0 b ) 


(Cl) 


The length f over which Yj ranges is specified as mR. If the effective noise- 
generating region is taken as greater than 5 diameters, then m > 10. Thus, for most 
values of 0 b and the second term is much less than the first term evaluated at 
Yj = L This disparity breaks down as Y^ — 0, but nevertheless the contribution of the 
second term in equation (Cl) to the integral in equation (30) will normally be minor, it is 
thought. The exceptions are 0 b » -0 and 0 & and 0 b near zero or 180°. In view of 
the mathematical complexity attendant on retaining the second term, this term will be 
dropped except near 0 b = -0 a . 2 The integral in equation (30) is then approximated as 

o )|cos0 b -cos0 a 

se 277 Jq 




= e -( 2 “f/ c o)| cos<, b- cosl) a| Y l (c2 

The remaining integral in equation (30) is then the integral of equation (C2) over Yj. If 
o>£ is treated as a constant 2 (^f), the integral is 


2 The subcase 0 b = -0 & of case n is precisely equivalent to subcase = 180° 

of case I. In this subcase only the second term of equation (Cl) is nonzero, and it leads 
to equations (46). Equations (47) are invalid for this subcase and are to be replaced by 
equations (46). 

2 The more realistic assumption GOf <=c Yj“l leads to an even simpler result. How- 
ever, the specification cuj = Constant is necessary for consistency with appendix B. 

(See last paragraph of appendix B for the reasons for the choice therein.) 
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*se “ 


1 - e 


-26 


26 


< w f) £ i 

6 = cos - cos 0 a 


A similar integral arises in the case of shear noise: 


(C3) 



1 - e 


-6 


6 


(C4) 


The ratios in equations (C3) and (C4) approach 0/0 as - 6 > a ; however, the limit 
is unity. This indeterminacy causes no difficulty with numerical evaluation if the case 
0^ = $ a is approximated by 0^ = 0 a + 0.001°. 

Equations (C3) and (C4) define the noncompactness factors arising from the integral 
in equation (30). For the complete reduction of equation (30) into equations (34), it is nec- 
essary to use the relation \j/= j0 b - 0j for the present case, where <p b = <p^. 
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REDUCTION OF THE SHEAR-NOISE INTEGRAL, EQUATION (42) 


The Yj - integration in equation (42) is approximated according to the ideas of 
appendix B (last paragraph) by replacing by an "effective" Yj -average ( w f)- 
Then the definition equation (44), allows the result to be rewritten more compactly as 


R sh(* a ’* b ’°) = ( p2 sh) max L S 0 [ c 


f~ 2 2 

COS 0 a COS 0^ 


1 -< w f>l A l 

+ cos 0 a cos sin sin 8 b sin^y - (p^sin yje 


dy (Dl) 


For case I: 0 & = = 8, cp a = 0, 

Rsh(5 a .x b ,0)>(p 2 sh L^.C’[ COs4e 


2 2 v • 1 "< w f>l A L 

+ cos 0 sin 9 sin^y - cp^jSin yje ay 


(D2) 


and 


A = — sin 6 sin(— | sin e 



which are equations (B4) and (B6), respectively. Upon making use of integral definitions 
introduced in appendix B, equation (D2) can be put in the form 



+ i cos 2 0 sin 2 0^cos cp^ - l) I^q sh ) 



(D3) 


where 


39 



APPENDIX D 


2/co f )R 

q sh = sin 9 sin 


(%) 
\ 2 ) 


‘NO'tlo 

r ( q S h) = i £ 


2 f 77/2 " q sh sin 6 


de 


ff /2 . 2 - q sh 


sin e 


sin e e 


de 


which are equations (B5), (B8), and (B9), respectively. 

Recall that the angle between x a and x^ (or between n a and n^) was defined 
as ip. Thus 

, a b a b a b 
cos + ^2^2 ^3^3 

= cos 0 a cos e b + sin 9 & sin 0 b (cos cp & cos <p b + sin <p a sin <? b ) 


= cos 0 a cos 6 h + sin 9 & sin cos^ b - cp^ 


(D4) 


For the present case (& a = © b = 9, <? a = 0), equation (D4) reduces to 


2 2 

cos ip = cos 9 + sin 8 cos cp-^ 


(D5) 


Equation (D3) may be reformulated in terms of i p as 


R s h(x a ,x b ,°)=(p 2 s hj max |{(' 


2 2 
+ cos 9 sin 9 


4 2 

[cos 9 + cos 9 cos 


*( q sh) 


■N 

j^( q sh) " I ( q sh) 


(D6) 
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SELF- SHEAR 


OTHER 

TERMS 


NOISE NOISE 


• THI S_P I CTURE I S GENERALIZATION OF RIBNER (1969) MODEL (ip = 0) FROM 

p l (1 MICROPHONE) TO p p u (2 MICROPHONES) 
d a b 

Figure 1.- Scenario for deriving two-microphone correlation p a p^. 
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OF TRAVEL TIMES FROM 
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A = f(0 AND OTHER ANGLES) 



SOURCE TERM INCLUDES 


dr 


a b 

4 u x u x “ 


“Ui 



• INSERTION OF x = A YIELDS CUSP- 

LIKE DECAY OF _ P a P b WITH V AND OTHER ANGLES 

• THIS IS AN EFFECT OF NONCOMPACTNESS OF THE SOURCE REGION 


Figure 2.- Effects of difference in travel times from source to p a and 

(noncompactness effect). 
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OF SOURCE REGION 


SOURCE-OBSERVER 
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Figure 3.- Physical basis of noncompactness effect (illustrated for vertical plane). 
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CASE I, WHEN \ 

0 = 0 = 9 °° 
a b / 


NONCOMPACTNESS 
OF SOURCE REGION 


SOURCE-OBSERVER 
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TRAVEL TIMES 



Figure 4.- p p h (i//) Source directivity x Noncompactness effect. 



Figure 5.- Geometry of problem for one microphone. 



Figure 6.- Geometry of microphone positioning. Case I: Both microphones lie on a 

cone coaxial with jet = 9; = 0, varies) . Case II: Both micro- 

phones lie on a meridian circle ( one at one at 0^) . 






U1 

o 



Figure 9.- Cartesian plots of normalized P a Pk (Vty for case I of figure 6 - comparison of 
experiment (Maestrello, ref. 7) with present theory for 6 = 90°. 




M. = 0.75 
EXPERIMENT 



Figure 10.- Cartesian plots of normalized P a P b (<P b , 0) for case I of figure 6 - comparison of 
experiment (Maestrello, ref. 7) with present theory for 6 = 120° and 135°. 








ANY PLANE PLANE p = e(t) PLANES Xj= CONST. 

<P= CONST. (HORIZONTAL) 


FIRST TERM SECOND TERM 

Figure 15.- Directional polars of the two components of instantaneous shear-noise 
pressure. The plane cp = e fluctuates randomly in orientation (fluctuating 
beam pattern concept). 
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